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Abstract 

We report astrometric observations of H2O masers around the red supergiant VY Canis Majoris (VY 
CMa) carried out with VLBI Exploration of Radio Astrometry (VERA). Based on astrometric monitoring 
for 13 months, we successfully measured a trigonometric parallax of 0.88 ± 0.08 mas, corresponding to a 
distance of 1.14 Iq'oq kpc. This is the most accurate distance to VY CMa and the first one based on an 
annual parallax measurement. The luminosity of VY CMa has been overestimated due to a previously 
accepted distance. With our result, we re-estimate the luminosity of VY CMa to be (3 ± 0.5) x 10 5 
Lq using the bolometric flux integrated over optical and IR wavelengths. This improved luminosity value 
makes location of VY CMa on the Hertzsprung-Russel (HR) diagram much closer to the theoretically 
allowable zone (i.e. the left side of the Hayashi track) than previous ones, though uncertainty in the 
effective temperature of the stellar surface still does not permit us to make a final conclusion. 

Key words: astrometry - masers (H2O) - stars: distances - stars: supergiants - stars: individual (VY 
Canis Majoris) - VERA 



1. Introduction 

Massive stars play an important role in the evolution of 
the Universe particularly in the late stages of their evo- 
lution. Through their strong stellar winds and supernova 
explosions, they inject mechanical energy into the inter- 
stellar medium (ISM) (Abbott 1982). Also, they are prin- 
cipal sources of heavy elements in the ISM. In spite of their 
importance, our understanding of late stages of massive 
stars' evolution is still poor. One of the reasons for this 
is that massive stars are extremely rare partly because 
of their short lifetime. Due to their small numbers, the 
properties of evolved massive stars are still uncertain. For 
instance, as discussed in Masscy (2003), Massey & Olsen 
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(2003) and Lcvcsque et al. (2005), there was a discrep- 
ancy between observed and theoretically predicted loca- 
tions of red supergiants, high-mass evolved stars, on the 
Hertzsprung-Russel (HR) diagram. Compared with stel- 
lar evolutionary models, observed red supergiants appear 
to be too cool and too luminous. 

VY Canis Majoris (VY CMa) is one of the most well- 
studied red supergiants. Similar to other red supergiants, 
the location of VY CMa on the HR diagram is also un- 
certain. For instance, Monnier et al. (1999), Smith et al. 
(2001) and Humphreys, Helton, & Jones (2007) obtained 
the luminosity of (2-5) x 10 5 L Q from the spectral en- 
ergy distribution (SED) assuming the distance of 1.5 kpc 
(Lada & Reid 1978) and effective temperature of 2800- 
3000 K based on the stellar spectral type (Le Sidaner & 
Le Bcrtre 1996). However, Massey et al. (2006) suggested 
that the above parameters would make VY CMa cooler 
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and more luminous than what current evolutionary mod- 
els allow and would place it in the "forbidden zone" of 
the HR diagram, which is on the right-hand side of the 
Hayashi track in the HR diagram (see figure 4). Massey 
et al. (2006) reexamined the effective temperature and ob- 
tained a new value of 3650 ± 25 K based on optical spec- 
trophotometry combined with a stellar atmosphere model 
and suggested that the luminosity of VY CMa is only 6.0 
x 10 4 L Q , which is probably a lower limit (Massey et al. 
2008). 

Levesque et al. (2005) determined the effective tempera- 
tures of 74 Galactic red supcrgiants based on optical spec- 
trophotometry and stellar atmosphere models. They ob- 
tained the effective temperatures of the red supergiants 
with a precision of 50 K. Their new effective tempera- 
tures are warmer than those in the literature. These new 
effective temperature values seem to be consistent with 
the theoretical stellar evolutionary tracks. However, the 
luminosities of the red supcrgiants still had large uncer- 
tainty due to possible errors in estimated distances. Since 
most of red supergiants are very far, it has been difficult 
to apply the most reliable trigonometric parallax method 
to their distance measurements. 

In the case of VY CMa, the currently accepted distance 
of 1.5 kpc (Lada & Reid 1978) is obtained by assuming 
that VY CMa is a member of the NGC 2362 star clus- 
ter and that the distance of NGC 2362 is same as that 
of VY CMa. The distance to the NGC 2362 was deter- 
mined from the color-magnitude diagram (Johnson et al. 
1961) with an accuracy of no better than 30 %. Since the 
luminosity depends on the square of the distance, the ac- 
curacy of luminosity estimation is worse than 50 %. The 
proper motions measured with H2O masers in previous 
studies (Richards. Yates, & Cohen 1998; Marvel et al. 
1998) have not been contradictory to the above distance 
value (1.5 kpc). When the H2O maser region around a 
star is modeled with a symmetrically expanding spherical 
shell, the distance to the star from the Sun can be inferred 
by assuming that the observed mean proper motion mul- 
tiplied by the distance should be equal to the observed 
mean radial velocity of the H2O masers. The proper mo- 
tion velocities in Richards, Yates, & Cohen (1998) are well 
consistent with the H 2 maser spectral line velocities at 
the distance of 1.5 kpc. Also, Marvel et al. (1998) derived 
the distance of VY CMa to be 1.4 ± 0.2 kpc. However, 
such a "statistical parallax" method depends on kinemat- 
ical models fitted to the observed relative proper motions. 
In fact, if we add rotation or anisotropic expansion to the 
simple spherical expansion model, the inferred distance 
value could vary beyond quoted error ranges. An accurate 
distance determination without an assumption is crucial 
for obtaining the right location of VY CMa on the HR 
diagram and for determining fundamental parameters of 
the star. 

Thanks to the recent progress in the VLBI technique, 
the distance measurements with trigonometric parallaxes 
have become possible even beyond 5 kpc (e.g., Honma et 
al. 2007). Since VY CMa has strong maser emission in 
its circumstellar envelope, we have conducted an astro- 



metric observations of H2O masers around VY CMa to 
measure an accurate parallax with VLBI Exploration of 
Radio Astrometry (VERA) and here we present the re- 
sults. 

2. Observations and Data Reduction 

We have observed H 2 masers (H2O 6i6-523 transition 
at the rest frequency of 22.235080 GHz) in the red super- 
giant VY CMa with VERA at 10 epochs over 13 months. 
The epochs are April 24, May 24, September 2, October 
30, November 27 in 2006, January 10, February 14, March 
26, April 21 and May 27 in 2007 (day of year 114, 144, 
245, 303, 331 in 2006, 010, 045, 085, 111 and 147 in 2007, 
respectively) . In each epoch, VY CMa and a position ref- 
erence source J0725-2640 [q(J2000) = 07h25m24.413135s, 
5(J2000) = -26d40'32. 67907" in VCS 5 catalog, (Kovalcv 
et al. 2007)] were observed simultaneously in dual-beam 
mode for about 7 hours. The separation angle between VY 
CMa and J0725-2640 is 1.059 degrees. The instrumental 
phase difference between the two beams was measured at 
each station during the observations based on the corre- 
lations of artificial noise sources (Kawaguchi et al. 2000; 
Honma et al. 2008). A bright continuum source (DA193 
at lst-3rd epochs and J0530+1330 at other epochs) was 
observed every 80 minutes for bandpass and delay calibra- 
tion at each beam. 

Left-handed circularly polarized signals were sampled 
with 2-bit quantization and filtered with the VERA digi- 
tal filter unit (Iguchi et al. 2005). The data were recorded 
onto magnetic tapes at a rate of 1024 Mbps, providing a 
total bandwidth of 256 MHz which consists of 16 x 16 
MHz IF channels. One IF channel was assigned to the 
H2O masers in VY CMa and the other 15 IF channels 
were assigned to J0725-2640, respectively. The correla- 
tion processing was carried out on the Mitaka FX correla- 
tor (Chikada et al. 1991). The spectral resolution for H2O 
maser lines is 15.625 kHz, corresponding to the velocity 
resolution of 0.21 km s _1 . 

All data reduction was conducted using the NRAO 
Astronomical Image Processing System (AIPS) package. 
The amplitude and the bandpass calibration for target 
source (VY CMa) and reference source (J0725-2640) were 
performed independently. The amplitude calibration of 
each antenna was performed using system temperatures 
measured during observation. The bandpass calibration 
was applied using a bright continuum source (DA193 
at lst-3rd epochs and J0530+1330 at other epochs). 
Doppler corrections were carried out to obtain the radial 
motions of the H2O masers relative to the Local Standard 
of Rest (LSR). 

Then, we calibrated the clock parameters using the 
residual delay of a bright continuum calibrator, which is 
also used for the bandpass calibration. The fringe fitting 
was made on the position reference source J0725-2640 
with integration time of 2 minutes and time interval of 
12 seconds to obtain residual delays, rates, and phases. 
These phase solutions were applied to the target source 
VY CMa and we also applied the dual-beam phase cali- 
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Fig. 1. The total power spectrum of the H2O masers in VY 
CMa obtained at Iriki station on April 24, 2006. 

bration (Honma et al. 2008) to correct the instrumental 
delay difference between the two beams. We calibrated 
the tropospheric zenith delay offset by using the GPS 
measurements of the tropospheric zenith delay made at 
each station (Honma, Tamura, & Reid 2008). After these 
calibrations, synthesized clean images were obtained. A 
typical synthesized beam size (FWHM) was 2 mas x 1 
mas with a position angle of -23°. We measured positions 
of the H2O maser features relative to the extragalactic 
source J0725-2640 with two-dimensional Gaussian fitting. 

3. Results 

The autocorrelation spectrum of the H2O masers 
around VY CMa is shown in figure 1 . The spectrum shows 
rich maser emission over LSR velocities ranging from -5 
to 35 km s _1 . Though there are variations in flux density, 
overall structures of maser spectrum are common to all 
the epochs, indicating that maser spots are surviving over 
the observing period of 13 months. 

To reveal the distribution of the H2O masers, we 
mapped the H2O maser features in VY CMa at the first 
epoch (April 24, 2006). We detected 55 maser features 
that have signal-to-noise ratio larger than 7 in two adja- 
cent channels. The distribution of the H2O maser features 
is shown in figure 2. The most red-shifted component and 
the most blue-shifted component lie near the apparent 
center of the distribution and the moderately red-shifted 
components on a NE-SW line of 150 mas. This distribu- 
tion well agrees with the previous observational results by 
Marvel (1996). 

Among the H2O maser features in VY CMa, one maser 
feature at a LSR velocity of about 0.55 km s _1 is ana- 
lyzed to detect an annual parallax. This feature, the most 
blue-shifted discrete component, is stably detected at all 
epochs. It is the third brightest maser component in the 
total-power spectrum, but in fact it is the second brightest 
component in the map. We note that while the brightest 
channel has multiple components, there is only one com- 
ponent in the channel at the LSR velocity of about 0.55 
km s _1 . Because of the simple structure as well as its 



strength, this velocity channel is most suitable for astro- 
metric measurements. Analyses for the other H2O maser 
features will be reported in a forthcoming paper. 

Figure 3 shows the position measurements of the 
0.55 km s _1 H2O maser component for 13 months. 
The position offsets are with respect to a(J2000.0) 
= 07 h 22 m 58 s .32906, <5(J2000.0) = -25°46'03".1410. 
Assuming that the movements of maser features are com- 
posed of a linear motion and the annual parallax, we ob- 
tained a proper motion and an annual parallax by the 
least-square analysis. The declination data have too large 
errors to be suitable for an annual parallax and a proper 
motion measurement. Therefore, we obtained the paral- 
lax of VY CMa to be 0.88 ± 0.08 mas, corresponding to 
a distance of 1.14 ~t ] 09 kpc using only the data in right 
ascension. This is the first distance measurement of VY 
CMa based on an annual parallax measurement with the 
highest precision. We estimated the positional uncertain- 
ties from the least-square analysis, and the values of the 
errors are 0.17 mas in right ascension and 0.68 mas in 
declination making reduced % 2 to be 1. 

We also determined the absolute proper motion in right 
ascension and declination. The absolute proper motion is 
-2.09 ± 0.16 mas yr _1 in right ascension and 1.02 ± 0.61 
mas yr _1 in declination. Compared with the proper mo- 
tion of VY CMa obtained with Hipparcos, 9.84 ± 3.26 mas 
yr _1 in right ascension and 0.75 ± 1.47 mas yr _1 in dec- 
lination (Perryman et al. 1997), there is a discrepancy in 
the proper motion in right ascension. Nearly 12 mas yr _1 
difference (more than 65 km s _1 at 1.14 kpc) seems too 
large even when we take into account possible relative mo- 
tion of the 0.55 km s _1 maser feature with respect to the 
star itself. The difference may be due to the fact that cir- 
cumstellar envelope of VY CMa has complex small-scale 
structure at optical wavelengths. This might seriously af- 
fect proper motion measurements with Hipparcos. 

The quantitative estimation of the individual error 
sources in the VLBI astrometry is difficult as previously 
mentioned in Hachisukaet al. (2006), Honma et al. (2007), 
and Hirota et al. (2007). Therefore, we estimated errors 
from the standard deviations of the least-square analysis 
to be 0.17 mas in right ascension and 0.68 mas in dec- 
lination. When we consider that the statistical errors of 
the position, 0.02-0.09 mas, are estimated from residuals 
of the Gaussian fitting, the large values of the standard 
deviations suggest that some systematic errors affect the 
result of astrometry. In the following, we discuss causes 
of these astrometric errors, although we cannot separate 
each factor of error. 

First, we consider the errors originating from the refer- 
ence source. The positional errors of the reference source 
J0725-2640 affect those of the target source. The position 
of J0725-2640 was determined with an accuracy of 0.34 
mas in right ascension and 0.94 mas in declination, re- 
spectively (Kovalev et al. 2007). Because these offsets are 
constant at all epochs, it dose not affect the parallax mea- 
surements. Also, when the reference source is not a point 
source, the positional errors of the target source could oc- 
cur due to the structure and its variation of the reference 
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Fig. 2. The distribution of the H2O masers in VY CMa obtained on April 24, 2006. Color, as shown on the scale on the right-hand 
side, represents the LSR velocities of the maser features. 
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Fig. 3. Results of the measured positions of the H2O maser spot at the LSR velocity of 0.55 km s 1 in VY CMa using 
J0725-2640 as a position reference source. The position offsets are with respect to a(J2000.0) = 07 h 22 m 58 s . 32906, <5(J2000.0) 
= -25°46'03" .1410. The left panel shows the movements of the maser spot in right ascension as a function of time (day of 
year). The right panel is the same as the left panel in declination. Solid lines represent the best fit model with an annual 
parallax and a linear proper motion for the maser spot. Dotted lines represent the linear proper motion (-2.09 ± 0.16 mas 
yr _1 in right ascension and 1.02 ± 0.61 mas yr _1 in declination) and points represent observed positions of maser spot with 
error bars indicating the positional uncertainties in systematic errors (0.17 mas in right ascension and 0.68 mas in declination). 
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source. However, since the reference source for our mea- 
surements is point-like and shows no structural variations 
between epochs, this is not likely to be the main source of 
the positional errors. 

Secondly, we consider the baseline errors originated 
from the positional errors of each VLBI station. The posi- 
tions of VERA stations are determined with an accuracy 
of 3 mm by the geodetic observations at 2 and 8 GHz every 
2 weeks. The positional errors derived from the baseline 
errors are 11 pas at a baseline of 1000 km with the base- 
line error of 3 mm. This error is much smaller than our 
astrometric errors. 

Thirdly, a variation of the structure of the maser fea- 
ture could be one of the error sources. Hirota et al. (2007), 
Imai et al. (2007), and Hirota et al. (2008) proposed the 
maser structure effect as main sources in their trigonomet- 
ric parallax measurements for nearby star forming regions 
Orion KL, IRAS 16293-2422 in p Oph East, and SVS 13 
in NGC 1333, correspondingly However, the effect docs 
not seem dominating in our case because, (1) the 0.55 km 
s _1 maser feature showed a stable structure in the closure 
phase, spectrum and map at all epochs of our observation, 
(2) this effect is inversely proportional to the distance of 
the target source and hence should be more than twice 
less significant for VY CMa as those in the above cases 
for a given size of the structure variation, (3) it is difficult 
to explain the large difference between astrometric errors 
in right ascension and declination by this effect. 

Finally, we have to consider the errors by the zenith de- 
lay residual due to tropospheric water vapor. These errors 
originate from the difference of path length through the 
atmosphere between the target and the reference sources, 
and generally larger in declination than in right ascen- 
sion. According to the result of the simulation in Honma, 
Tamura, & Reid (2008), the positional error by the tro- 
pospheric delay is 678 pas in declination when the atmo- 
spheric zenith residual is 3 cm, the declination is -30 de- 
grees, separation angle (SA) is 1 degree, and the position 
angle (PA) is degrees. This is well consistent with our 
measurements. Therefore, the atmospheric zenith delay 
residual is likely to be the major source of the astrometric 
errors. 

4. The Location on the HR Diagram 

We successfully detected a trigonometric parallax of 
0.88 ± 0.08 mas, corresponding to a distance of 1.14 Ig'Jg 
kpc to VY CMa. Compared with the previously accepted 
distance 1.5 kpc (Lada & Reid 1978), the distance to VY 
CMa became 76 %. Since the luminosity depends on the 
square of the distance, the luminosity should become 58 
% of previous estimates. Hence, here we re-estimate the 
luminosity of VY CMa with the most accurate distance. 
The luminosity can be estimated as follows: 

L = 4nd 2 F bol , (1) 

where L is luminosity, d is distance and F^oi is the bolo- 
metric flux. To obtain Fb i, we used the SED of VY CMa. 
The data are based on HST optical images and near-IR 
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ground based images in Smith et al. (2001) and IRAS 
fluxes from 25 to 100 pm. The Ftoi is obtained by inte- 
grating the observed fluxes. The estimated luminosity of 
VY CMa with our distance is (3 ± 0.5) x 10 5 L . 

We re-estimate the luminosities of VY CMa in the pre- 
vious studies with our distance. Le Sidaner & Le Bertre 
(1996) obtained a luminosity of VY CMa to be 9 x 10 5 
L Q from the SED at a kinematic distance of 2.1 kpc. With 
the distance of 1.14 kpc, the luminosity of Le Sidaner & Le 
Bertre (1996) became 2.6 x 10 5 L . Smith et al. (2001) 
also estimated a luminosity of VY CMa from the SED and 
their luminosity is 5 x 10 5 L Q at a distance of 1.5 kpc, 
which is is revised to be 3.0 x 10 5 L Q using our distance. 
These luminosities are well consistent with each other. 
Massey et al. (2006) estimated an effective temperature 
of VY CMa to be 3650 K, which fitted the MARCS at- 
mosphere model to the observed spectrophotometric data. 
They estimated an absolute V magnitude My using the 
observed V magnitude, currently known distance of 1.5 
kpc and a visible extinction Ay- Since the MARCS model 
presents bolometric corrections as a function of effective 
temperature, they obtained the luminosity of VY CMa 
to be 6.0 x 10 4 Lq, much lower than our value derived 
above. When they adopt our distance, their luminosity 
would be increased. The accurate distance is essential to 
estimate the accurate luminosity. 

To place VY CMa on the HR diagram, we adopted an 
effective temperature from the literatures. As we already 
mentioned, Massey et al. (2006) estimated an effective 
temperature of VY CMa to be 3650 K based on the ob- 
served spectrophotometric data. With the previously ac- 
cepted spectral type of M4-M5, Le Sidaner & Le Bertre 
(1996) obtained an effective temperature of 2800 K and 
Smith et al. (2001) also adopted an effective temperature 
of 3000 K. Since the effective temperature does not depend 
on the distance, our measurements cannot judge which ef- 
fective temperature is correct. 

Although we cannot estimate the effective temperature 
of VY CMa with our measurements, when we adopt the 
effective temperature of 3650 K from the MARCS atmo- 
sphere model (Massey et al. 2006), the location of VY 
CMa on the HR diagram is determined as the filled square 
in figure 4. For comparison, we also show VY CMa's loca- 
tions on the HR diagram obtained in the previous studies 
as filled circles. Our results suggest that the location of 
VY CMa on the HR diagram is now consistent with the 
evolutionary track of an evolved star with an initial mass 
of 25 Mq. Also, re-scaled luminosity values of Le Sidaner 
& Le Bertre (1996) and Smith et al. (2001) imply much 
closer locations to the 25 M Q track than their original 
ones which were deeply inside the "forbidden zone". On 
the other hand, the lower limit of luminosity of 6.0 x 10 4 
L Q in Massey et al. (2006) would be consistent with 15 M Q 
initial mass of VY CMa. According to Hirschi, Meynet, 
& Maeder (2004), there is an order of magnitude differ- 
ence in lifetimes between 15 M Q and 25 M© in the initial 
mass. The improvement in mass and age values should af- 
fect statistical studies on evolution of massive stars such 
as the initial mass function (Salpeter 1955). 
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On the other hand, there is an argument against the 
effective temperature in Massey et al. (2006). Humphreys 
(2006) suggested that the spectrum in Massey et al. (2006) 
are more like their M4-type reference spectrum than M2- 
type reference spectrum. Humphreys (2006) also pointed 
out that the modeling applicable for "standard" stars 
without mass-loss (MARCS) is simply not valid for an 
object like VY CMa. When we adopt the previously 
accepted spectral type, instead of M2.5I (Massey et al. 
2006), the effective temperature is 3000 K (Smith et al. 
2001). The location of VY CMa on the HR diagram is 
shown as the open square in figure 4. In this case, the po- 
sition on the HR diagram is still not consistent with the 
theoretical evolutionary track. Table 1 summarized our 
adopted parameters for VY CMa. 

The accurate distance measurements of red supergiants, 
which provide true luminosities, will greatly contribute 
to the understanding of massive star evolution, though 
there is still uncertainty in the effective temperature on 
the stellar surface. 

5. Conclusion 

We have observed the H2O masers around the red su- 
pcrgiant VY CMa with VERA during 10 epochs spread 
over 13 months. Simultaneous observations for both H2O 
masers around VY CMa and the position reference source 
J0725-2640 were carried out. We measured a trigono- 
metric parallax of 0.88 ± 0.08 mas, corresponding to a 
distance of 1.14 kpc from the Sun. It is the first 

result that the distance of VY CMa is determined with 
an annual parallax measurement. There had been over- 
estimation of the luminosities in previous studies due to 
the previously accepted distance. Using the most accu- 
rate distance based on the trigonometric parallax mea- 
surements and the bolometric flux from the observed SED, 
we estimated the luminosity of VY CMa to be (3 ± 0.5) 
x 10 5 L Q . The accurate distance measurements provided 
the improved luminosity. The location of VY CMa on 
the HR diagram became much close to the theoretically 
allowable region, though there is still uncertainty in the 
effective temperature. 

We are grateful to the referee, Dr. Philip Massey, for 
helpful comments on the manuscript. We would like to 
thank Prof. Dr. Karl M. Menten for his invaluable com- 
ments and for his help improving the manuscript. The 
authors also would like to thank all the supporting staffs 
at Mizusawa VERA observatory for their assistance in ob- 
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Table 1. Adopted parameters for VY CMa 



Parameter 


Value 


Note 


Parallax 


0.88 ± 0.08 mas 




Distance 


1.14 i° H kpc 




Luminosity 


(3 ± 0.5) x 10 s L Q 




Mass 


25 Mq 


Meynet & Maeder (2003) 


Temperature 


3650 ± 25 K 


Massey et al. (2006) 




3000 K 


Smith ct al. (2001) 




Fig. 4. The various locations of VY CMa on the HR diagram. The filled and the open squares represent our results. The luminosity 
of our result is calculated using the distance based on the trigonometric parallax measurements and the bolometric flux from the 
SED. The filled square adopted the effective temperature of 3650 K (Massey et al. 2006) and the open square adopted 3000 K(Smith 
ct al. 200f). The circles "1," "2" and "3" represent the results of Le Sidaner & Le Bertre (1996), Smith et al. (2001) and Massey et 
al. (2006), respectively. The evolutionary tracks are from Meynet & Maeder (2003). 



